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Abstract 

A celestial body produces a gravitational moment about the mass center of a 
small orbiting body, which affects the orientation of the smaller body. Each 
zonal harmonic in the gravitational potential of a celestial body is shown to 
make a contribution to the gravitational moment which can be expressed in a 
recursive vector-dyadic form. A formal derivation is presented, followed by an 
example in which the result is employed in obtaining the contribution of the 
zonal harmonic of 2nd degree. The contribution of the zonal harmonic of 3rd 
degree is also reported. 

Introduction 

The gravitational moment about the mass center of a body in orbit about a 
celestial body has an important effect on the orientation of the orbiting body. 
The more misshapen the celestial body, and the less uniform its mass distri- 
bution, the more involved is the calculation of the gravitational moment (and 
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force) it exerts. Situations in which it might be important to calculate accu- 
rately the gravitational moment include the design of spacecraft for expeditions 
to asteroids, comets, and the moons of Mars. 

In Ref [1], a method for obtaining a vector-dyadic expression for the moment 
exerted about a small body’s mass center by an oblate spheroid was set forth. 
The derivation of that expression made use of a gravitational potential written 
in terms of the zonal harmonic of 2nd degree. When gravitational potentials 
containing zonal harmonics of degree 2 or greater are considered, each zonal 
harmonic makes a contribution to the gravitational moment. 

Recorded below is a vector-dyadic expression for the contribution of a zonal 
harmonic of degree n to the gravitational moment, produced by a body, about 
the mass center of a small body. As is the case with all vector-dyadic expres- 
sions, this result is basis independent, — that is, the vectors and dyadics can be 
expressed in any convenient vector basis. 

The equation given below is recursive: the contribution to the gravitational 
moment from the zonal harmonic of degree n is a function of the moment con- 
tributions from the zonal harmonics of degree n — 1 and n — 2. The equation 
contains Legendre polynomials and derivatives of Legendre polynomials of de- 
gree n — 1 and n — 2. The Legendre polynomials, as well as their derivatives, 
can, themselves, be generated by means of recursion formulae. 

As an example, the contribution to the gravitational moment from the zonal 
harmonic of degree 2 is worked out. The contribution of the zonal harmonic of 
degree 3 is also given. 

Assertion 

Figure 1 shows a small body B in the presence of an axisymmetric body E. 
The distance between £*, the mass center of B, and E* , the mass center of 
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E, is assumed to exceed the greatest distance from jB* to any point of B . The 
system of gravitational forces exerted by E on B produces a moment M about 
£*, and M can be written as 

o OO 

M= |rx /.r + £M„ (1) 

n=2 

where M n is the contribution of the zonal harmonic of degree n and can be 
obtained by using the recursion relation 


Mn = [ (2n - 1} + 2)Pn - 1 + 3( * • ~ (2n - 2)p "- 2 1 


(to x / • r + # x Jtn) 


+ [(4n - 4) [(n + 1 )P n -i + (r • n)P^_ 2 ] 

— (2n — l)(r • n) [(4n + 8)P n _i 4- 4(f • n)P'_ 1 ]Jf x J • f 

— j^2(2n — J x L • 

■ 2u — 1 J n Re / ~ * \ n/r n — 1 J n f Re \ „ ^ 

+ — 3^7-R (r,n)M ”-‘-— J^Vr) m “- 2 


where; \i is the gravitational parameter of E , J n is the zonal harmonic coefficient 
of degree n, Re is the mean equatorial radius of E, R is the distance from E* 
to B*, P n is the Legendre polynomial of degree n and argument Sa, Sa is the 
Sine of A, the latitude of B* (Sa = r *to), P. ^ is the first derivative, with respect 
to its argument, of P n , r is the unit position vector from E * to B* , to is the 
unit vector in the direction of the axis of symmetry of £*, and / is the inertia 
dyadic of B relative to B*. 
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Figure 1: Body B in the Presence of Body E 

Derivation 

The gravitational forces exerted by a body E on a small body produce a moment 
M about the mass center £* of B. M is given approximately by equation 
(2.18.1) of Ref [2], 

M ==-/. VW(R) (3) 

where R is the position vector from E* to B* and V denotes differentiation with 
respect to the vector R. Section 2.9 of Ref [2] contains a thorough explanation 
of how one differentiates with respect to a vector. The definition of the cross-dot 
product, •, appears on p. 156 of Ref [2]. The gravitational potential, of E is 
symbolized by V. 

Equation (2.13.14) in Ref [2] deals with the gravitational potential of an 
axisymmetric body and contains an infinite series of zonal harmonics. For a 
particle of unit mass coincident with 
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(4) 


V 


R 


1 - £ 

n—2 v ' 


where Sa, the argument of P n , is equal to the sine of the geographic latitude of 
B *. Eq. (4) can be simplified to 


V 


R 




(5) 


n= 2 


when the contribution to the gravitational potential of the zonal harmonic of 
degree n, V n , is defined as 


A fiRE n j p 

RP+ 1 nn 


( 6 ) 


Ref [1] shows that — J • W(fi/R) = (3n/R 3 )r x J*r. Hence, eq. (3) can 
be rewritten as eq. (1), 

o 00 

M= -£r x/.f + J^M n (7) 

n~ 2 

so long as the contribution to the gravitational moment of the zonal harmonic 
of degree n, M n , is defined as 


M n = -I 2 VVV n 


( 8 ) 


The Legendre polynomial of degree n, P n (z), is expressed recursively in 
equation (8.71) of Ref [3] in terms of Legendre polynomials of degree n — 1, 
n — 2, and their argument, x, for n > 2, as follows: 


P n (x) = —[(2n - l)xP n -i(x) - (n - l)P„_ 2 (a:)] (9) 

71 

Eq. (9) can also be produced with m = 0 in formula I of Table 1 in Ref [4]. Sub- 
stituting from this recursion relation for the Legendre polynomials into eq. (6) 
leads to a recursive expression for V n : 


Vn 


2n — 1 J n Re 

TX Jn — 1 R 1 


S a K„_i- 


n - 1 Jn f Re \ 
Tl Jn—2 V R ) 


2 

V n -2 


( 10 ) 


421 



Eq. (8) requires that a dyadic be formed by differentiating V n twice with 
respect to R. The first derivative of V n with respect to R yields a recursion 
relation for the vector VV n , the contribution of the zonal harmonic of degree 
n to the gravitational force exerted by E on a particle of unit mass coincident 
with B* . 


VF n 


2n — 1 J n 


■Re 


n J n - 1 

»- 1 J n 


Vn- 

R 


ivS A + %W n _t - 


R 


S aK,-1 

R 3 


R 


Tl Jn — '. 


■ R b 2 


2 f Wn— 2 

R? 


W n -2 ] 

~'rT R \ 


(U) 


The second derivative of V n with respect to R yields a recursion relation for the 
symmetric dyadic VVF n> 


VW n = 


+ 


2n — 1 


n 


n— 1 ^ K it 


Jn- 


VF„-iVSa , Sa.,^ S\W n -iR 

R + R n ~ 1 R? 

Vn-iRV Sa Sa RVVn-! 3 SaK „-1 RR 

R 3 R 3 + R 3 

n — 1 J n D 2 f ^^Kj -2 2 VVJj_ 2 -R 

n J n - 2 i. R? Rj 

2RVV n ^ 8V n _ 2 RR 2V„_ 2rr ' l 
RJ + R G R4 ~i 


V^VSa-R 

R 3 


SxK- 

R 3 



( 12 ) 


where {/ is the unit dyadic. 
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Cross-Dot Products 


We now perform the cross-dot product with dyadics / and VW n , making use 

X 

of the right hand side of eq. (12) and the cross-dot identity I « U_ .= 0, which is 
set forth in eq. (19) of Ref [1], obtaining 

L • vw„ = ~ ~ 1 ?Z- R E {¥*j±l S VVSa - • (VSaR+KVSa) 

+ • VVFn-i + i/ • (VSxVK,-! + V^-iVSa) 

- » (VK-ii? + ■RVV'n-i) + — —■ I • RR - o} 

- — J 5 VVKn-2 + • RR— 0 

72 Jn—2 1 R I 2 It 6 

- J ? I*(VV n _ 2 H + BVy„_ 2 )} (13) 


In order to carry out the cross-dot products with I and the other dyadics on 
the right side of eq. (13), we will express these dyadics in terms of R and n, a 
unit vector parallel to the axis of symmetry of E . 

The sine of A can be expressed as sin A = (R* n)/R = r • n so that the first 
derivative of S\ with respect to R is 


Rmn n (Rmri)R 
~ R~~ = R " 


(14) 


and the second derivative of Sa with respect to R is 

VVSa = 3(R ^ )Jlfi - -p [nR + Rh + (R. n)U] (15) 

Eq. (20) of Ref [1] is a derivation of a cross-dot identity which will be used 
repeatedly throughout the sequel: For any dyad uv composed of vectors u and 
u, it is shown that / • uv = —u x / • v By making use of this identity, and 
eq. (15) above, one can write 

I • VVSa = i (n x / • f + r x / • n) + 0 — x I • r (16) 

l? 2 “ 
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where r is a unit vector in the direction of R. By recalling the definition in 
eq. (6), one can evaluate the first cross-dot product on the right side of eq. (13). 

# VVS A == 

—lxJn-iRE n ~ l x L* t + v x I.#n — 3(r *n)r x J«f] (17) 

The dyadic required for the second cross-dot product on the right side of 
eq. (13) is easily constructed by using eq. (14): 

VS \R + RVS\ =nf+ rfi — 2{r • n)r r (18) 

Thus, 

Yjz±li(VS x R+RVS x ) = 

—fiJ n -iRE n ~ 1 ~^ [2(r • h)r x !• f — n x /• r — ** x £*n] (19) 

Replacing n with n — 1 in eq. (8) allows the third cross-dot product on the 
right side of eq. (13) to be immediately expressed in terms of M n _i . 

• VW„_! = -^M,: (20) 

The fourth cross-dot product on the right side of eq. (13) contains the quan- 
tity W n _i. Temporarily allow x to be the argument of P n -i [See eq. (6)] and 
write the derivative of K*_i with respect to R as 

[^r-^ Pn " l(a:)Va: " ! %^ p ] < 21 > 

The first derivative of a Legendre polynomial P n with respect to its argument 
is often denoted by P^. For n > 2, a useful recurrence formula for P„ can be 

found in problem 8-9 of Ref [3], p. 393, or in formula I of Table 1 in Ref [4] 

(with m = 1): 

K{X) = 7^1 f(2 ” “ l > P n~i( X ) - nP «-2(*)} ( 22 ) 


424 



Since the argument of P n -i is known to be Sa, we can make use of eq. (14) to 
rewrite eq. (21) as 


T 7 T/ - _ ..r D n -i[Ci fn (R*n)R] nP„_i „\ 

(-^[5 S5— \--^» R ) 


(23) 


The sum of two dyadics, formed by juxtaposing the vectors VSa and VV r n _ 1 in 
opposite order, yields the symmetric dyadic 


VSaVV„., + W„_,VSa = 


Consequently, the fourth cross-dot product on the right side of eq. (13) can be 
expressed as 


n— 1 


xl • (VS A VF^! + VK-iVSa) ^ -iiJn-iRz 

{ R™+ 3 ‘ # x/»n+nx/*f) - n x 7# n - (r • Ti) 2 ^ x 7 • rj 

TlP 1 

— [2(r • n)f X 7 • r — n x 7 • f — r x J«n]| (25) 

The dyadic required for the fifth cross-dot product on the right side of 
eq. (13) can be constructed rather easily by employing eq. (23), which yields 

-{i 3 n -\RE n ^ X | +'** “ 2(f • n)rr] - 2n ^ z - rr j (26) 


so that 


• (W^R + RVV n -!) = 


/iJn-i-Rf;" 1 ^ f ftn+ 3 1 [2(* • w)f x / • r — n x 7 • r — r x / • n] 
xl»r| (27) 


2n(f • n)P„_ i . 
+ =r— rs r 


J ?«+ 3 


425 



The sixth cross-dot multiplication which must be performed in order to 
x 

obtain I • WV n is one of the easiest to carry out. That is, 

3S ^ n -I »RR= ZnJn-iRE n ~\r . n)~^T xl.r (28) 

The seventh cross-dot product on the right side of eq. (13) can be expressed 
in terms of M n ~ %. Replacing n with n — 2 in eq. (8), we get 

T/ • VWn-2 = -pMn-S (29) 

The eighth cross-dot product to be evaluated is simply 

L • RR = ^Jn- 2 RB n ~ 2 ^r x 7 . f (30) 

The dyadic required for the final cross-dot product is similar to that needed 
for the fifth cross-dot product, 


VV n ^ 2 R + RVV n „ 2 = 

-liJ n - 2 RE n ~ 2 j^rf [fir + rn - 2 (f • n)r f] - ^ - R ^ n ~ 2 rr |(31) 
so that 


Jjl* (VV r „_ 2 fl+RVV r „_ 2 ) = 


f 2P f 

—flJn-2RE n ~ 2 S ftn +3 [2(f e h)r X / • £ — 72 X / • £ — £ x / ® ra] 


4(n-l)P„_ 2 

+ R n+3 


l • f 


(32) 


Substituting from eqs. (17), (19), (20), (25), (27) - (30), and (32) into (13) 
and then into (8) leads to eq. (2), which is a recursive vector-dyadic expression 
for the contribution of the nth zonal harmonic to the gravitational moment. 


Examples 

In order to demonstrate the use of eq. (2), we will use it to obtain M 2 , the 
contribution to the gravitational moment from the zonal harmonic of degree 2. 
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The two required Legendre polynomials are Pi(Sa) = S\ = r • ri and 
Pq( S\) == 1 . Legendre polynomials of degree greater than or equal to 2 can 
be obtained recursively by using eq. (9). The two derivatives of Legendre poly- 
nomials which will be needed are P[( Sa) = 1 and Pq(Sa) — 0. Derivatives with 
respect to the argument of Legendre polynomials can be generated with the 
recursion formula (22) for n > 2. 

Eq, (2) also requires knowledge of Af 0 and in order to produce M 2 . 
Eq. (8) is helpful in developing expressions for Mq and M\. 

The Legendre polynomial of degree zero is equal to 1, regardless of its argu- 
ment, and the scalar Vq [See eq. (6)] is 

Vo = -fiJoPo/R = —fiJo/R (33) 


M>, 


The dyadic formed with Vq is then 

VV V 0 = VV(-/iJo/i?) = V(nJ 0 R/R : 3 ) = ^(u -3 rr) (34) 
Eq. (8) tells us that 


Mo = x /• r 


(35) 


Jo is an undefined constant, but the coefficient of Mo in eq. (2) contains J 0 in 
the denominator. Hence, a numerical value of J 0 is not required for constructing 
M 2 . 

A similar process leads to Mi. The value of the Legendre polynomial of 
degree 1 is identical to the argument, so the scalar Vi is 


t/ T r. p t T D (R*n) 

Vi = -nJiRE-jg = -uJiRe - r 3 


(36) 


The dyadic formed with V) is 

VVVi = -hJiR e w ( R *^ = -hJxRbV 


h 3(R*n) 


]= 

-IxJiRe | 15( ~| y — RR -j^[Rn + hR + (R* n)&]} (37) 
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so that 


M i = ***« ■ [rix/®f+r X I ® n — 5(f ® fi)r x 7 • f] (38) 

Like Jo, the constant J\ is undefined, and unneeded for the purpose of obtaining 
M 2. Note that Mo and M\ do not represent contributions to the gravitational 
moment, but are required to begin the process of recursion which will generate 
moment contributions beginning with M 2 . 

By substituting from eqs (35) and (38) into (2), we arrive at the following 
result with n = 2: 

M 2 = ®n)](n x /»r + f x 7®n) 

4- |l5 — 105(r • n) 2 j r x I • r — 6n x J • n | (39) 

If eq. (1) of Ref [1] is expressed as M = (3 ji/R 3 )r x /•£ +.M2, it can be seen 
that M2 from Ref [1] is identical to eq. (39), above. 

The contribution M 3 can be obtained in a similar manner, using the values 
of P 2 (S A ), JY( S A ), PR$x), P[(S a ), M 2 , and M,. 

M$ == ^^^6 { [^15(r • ft) 2 — 45 j (nx/if-ff x J • n) 

+ J;315(f • n) — 945(r • n) 3 r x /• r — 90(f • n)n x / • n| (40) 

Conclusions 

A recursive vector-dyadic expression for the contribution of a zonal harmonic 
of degree n to the gravitational moment about the mass center of a small body 
can be obtained by a procedure which involves differentiating a celestial body's 
gravitational potential twice with respect to a vector. The recursive property 
of the result is a consequence of taking advantage of a recursion relation for 
Legendre polynomials that appear in the gravitational potential. When a celes- 
tial body's gravitational potential includes zonal harmonics, the vector-dyadic 
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expression above is useful for calculating their contributions to the gravitational 
moment. The contribution of the zonal harmonic of degree 2 is consistent with 
the gravitational moment exerted by an oblate spheroid. 
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